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We study numerically a disordered model that interpolates among the Sherrington-Kirkpatrick 
mean field model and the three dimensional Edwards- Anderson spin glass. We find that averages 
over the disorder of powers of the overlap and of the full P{q) are smooth, and do not show any 
discontinuity. Different lattice sizes are used to provide evidence for a smooth behavior of disorder 
averages in the thermodynamic limit. Quantities defined on a given realization of the disorder show 
a chaotic behavior. Our results support the validity of a Replica Symmetry Breaking description of 
finite dimensional models. 
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The relation among the mean field spin glass model 1 1 ] 
(SK) and its finite dimensional, realistic counterpart |2] 
(EA), is a debated and interesting subject. If on one side 
the Parisi replica symmetry breaking (RSB) solution [^[^] 
of the mean field theory is getting a stronger acceptance 
even on a rigorous ground on the other side the 

relation of RSB theory with finite dimensional models is 
less clear. 

Numerical simulations (sec Q and references therein) 
seem to confirm the validity of the Parisi picture in 3 
and 4 spatial dimensions D, but the droplet model ap- 
proach suggests a different behavior of finite dimen- 
sional models. Recent analytical, rigorous and heuristic 
work can be seen as hinting for potential troubles in the 
application of RSB to finite dimensional models |10|, or 
as confirming the suitability of such a description |ll,|l2|. 

If the droplet model description ^ applies to finite 
dimensional spin glasses the mean field model (that we 
assume as described by the Parisi RSB Ansatz) has to be 
a singular case: at variance with the usual Ising model 
the case of = oo would be different in nature from 
the case of any finite dimensionality. All the exceptional 
features that appear in the description of the mean field 
model solution would disappear for any large but finite 
dimension D. Typical equilibrium configurations of the 
SK and the EA models would have to be very different. 
On the contrary validity of the RSB ansatz in the finite 
dimensional case would imply that phenomena like the 
existence of a phase transition in non-zero magnetic field, 
and some form of ultrametric organization of the states 
appear also for finite dimensional models. 

Here we analyze numerically this issue, and show that 
there are no discontinuities in the typical behavior of the 
system (clearly the behavior of a single system turns to 
be, as expected, very discontinuous, while average over 
the disorder, i.e. the quenched state, does not show any 
abrupt discontinuity). There are many contexts where 
this result can be relevant. For example Guerra |^ has 
shown rigorously the validity of relations like Q 



'-E{qU) + -E{ql,) 



and 



E{ql2qlA) = l,E{ql^) + -E{ql, 



(1) 
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where q^^b if the overlap ^t'^'l among replica a and 

replica h of the system, and E{-) denotes the expecta- 
tion value over the thermal noise and over the quenched 
disorder. The proof holds for mean field, and, under a 
continuity hypothesis, for finite dimensional models. Nu- 
merically one finds ||l2| that in iD relations like (0J|) 
are satisfied at better than one part over thousand (with 
E{q^) different from zero in the spin glass phase). This 
is one of the reasons for which checking the continuity of 
the limit of finite dimensionality is an important task. 

Another relevant example where a weak perturbation 
is used to derive important relations among equilibrium 
observables and dynamic quantities that can be measured 
in experiments is reference [|l^. Here we will analyze a 
system where the perturbation does not need to be small: 
we want to understand the nature of the relation among 
the mean field model SK and the finite dimensional EA 
Ising spin glass. In our case starting from a mean field 
interaction involving with p > 2 spin would not give a 
smooth limit to a two spin interaction three dimensional 
model, since the mean field starting point would be a one 
step RSB P{q). 

We have defined and studied a model that interpolates 
from mean field to 3D. The Hamiltonian of our model is 
a function of the parameter e, which ranges from zero to 
one: 
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All the couplings Jij are binary, and take with prob- 
ability one half the values ±1. The three dimensional 
structure of the system is given by the term 



H: 



3D 
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where the sum runs over first neighboring sites on a sim- 
ple cubic lattice in three dimensions, with N = sites. 
The mean-field, Sherrington Kirkpatrick like contribu- 
tion is given by the term 



(5) 



where the sum runs over all site couples of the lat- 
tice. The quenched coupHngs {J^^\ J^^^} are indepen- 
dent. For e = we recover the ZD EA spin glass, while 
for e = 1 we recover the SK mean field model. 

We have focused our analysis on the probability distri- 
bution of the overlap, P'^{q), as a function of e, averaged 
over the quenched disorder, and on the P}{q) for fixed 
realization of the disorder { J'"'^-', J*^^^}. We have also an- 
alyzed in detail the second, fourth and sixth momentum 
q^, q'^ and q^ . 

As we will discuss our results clearly show that P^{q) 
is a smooth function of e in all the range of e ranging 
from zero to one, for, say T ~ 0.7Tc. On the contrary 
for a given disorder sample the shape of P}{q) depends 
dramatically on e. The expectation values of the first 
even powers of the overlap all have a smooth dependence 
over e. Larger lattices show a behavior consistent with 
this picture, and we do not see any singular behavior 
developing. In the clear limits of a numerical simulation 
we can state are observing a smooth transition from the 
mean field to the finite dimensional model. 

The Parallel Tempering Monte Carlo approach (for a 
review see for example has been crucial for ther- 

malizing our systems down to T ~ 0.7Tc. Since Tc = 1 
for the mean field model and Tc ~ 1.11 [|l^ is a reason- 
able estimate for the 31? system with binary couplings we 
have analyzed our data as a function of e at fixed T — 0.7. 
Renormalizing the T value as a function of e could have 
been in principle a more exact procedure: we have de- 
cided not to do so not to add a further degree of freedom 
(the scaling of the temperature) that would have in it- 
self potentially added a chaotic behavior, without adding 
much to the interpretation of the result. We are looking 
for smoothness of the observables (on the state averaged 
over the disorder) and not for precise extrapolations, so 
we can afford ignoring the small renormalization of Tc 
(that we could only control approximately in any case). 
All the data we will present in the following will be at 
T = 0.7. 

This numerical simulation is expensive, since for all 
e values we pay the price {N'^) of a mean field simula- 
tion (this is only the scaling of the number of operations 
needed, and does not include an estimate for the incre- 
ment of the needed computer time due to critical slowing 
down, so that it has to be read as a lower bound to the 
increase of the computer time): the time taken from a 
lattice sweep scales like the volume squared instead than 
like the volume. 

At L = 8 we have been able to analyze 150 disorder 
samples, each with lie values (from to 1.0 with steps of 



0.1), while for L = 12 we have 11 disorder samples with 
the same e values. As we have said this is not a small 
computational effort because of the factor volume added 
by the mean field interaction term and the factor eleven 
for the different e values that multiply the total computer 
time. We have used a multispin coded program (for the 
long range part of the interaction). The simulations have 
taken a few years of a nowaday typical workstation. 

The Parallel Tempering Monte Carlo updating scheme 
worked effectively at all our e values for both L = 8 and 
i = 12. For L = 8 we have used 13 copies of the sys- 
tem, with temperatures starting with 0.7 and increasing 
with steps of 0.05 up to 1.3, while at i = 12 we have 
used 26 copies with the same minimum value of T, up 
to T = 1.95 with the same step. We ran 100.000 lattice 
sweeps (and tempering updates) for each disorder real- 
ization and e value both at L = 8 and at i = 12. At all 
values of the parameters we have discarded the first half 
of the Monte Carlo sweeps for sake of thermalization. At 
each new e value we restarted the system from a random 
configuration, even when running the same disorder sam- 
ple (in order to avoid any bias toward continuity). All 
thermalization tests |lj] were satisfactory: the T swap 
acceptance ratio was included among 0.3 and 0.7, and 
each copy of the system had eventually spent time in all 
allowed T values. 

We have used the two different lattice sizes {L — 8 and 
L = 12) to get control over the finite size scaling of our 
results. It turns out that all the features we have revealed 
and we will discuss now are stable when increasing the 
lattice size, and make us confident we are observing an 
asymptotic behavior. The L = 12 results have a larger 
error than the L — 8 ones (since in this case the N'^ scal- 
ing of the computer time needed for the simulation has 
limited us to 11 samples), but they confirm the behavior 
we are observing at L = 8. 

In figure (^ we show the average probability distribu- 
tion of the overlap q, P{q), for the L = 8 lattice, as a 
function of e (and q), at T = 0.7. The e axis goes from 
zero to one, the overlap axis from —1 to 1. The behavior 
of the function P'^{q) as a function of e is very smooth. 
It is important to notice that the value of P{q = 0) at 
e = is largely different from zero (in a definitely signifi- 
cant way: this is very well known, see for example Q and 
references therein) . Continuity of P'^ (q) as a function of e 
seems very clear from figure (^. For L = 12 we get a very 
similar picture (see also the value of g^(e) that we will 
show later on), and things stay as smooth as they are in 
figure (^). The non trivial probability distribution that 
one finds in the mean field approach survives smoothly 
down to the finite dimensional, realistic system. 

We do not expect continuity for a given realization of 
the quenched disorder: the mean field picture would sug- 
gest that the typical equilibrium configurations are highly 
unstable under variations of e Q|. Our results confirm 
this point of view in a very nice way. We show in fig- 
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and (t 



) two Pj{q) for two typical realizations of 



the quenched disorder. In both figures one sees that the 
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shape of P){q) varies completely with e. In figure ^ the 
mean field model is in a configuration that carries all the 
weight close to g = 0. Already at e = 0.9 the system is 
completely different, and Pj'^{q) has two maxima. No- 
tice that the level of symmetry of these Pj (q) is a good 
check of the level of thcrmalization of our Monte Carlo 
tempered runs. The pure SD system has in this case a 
broad distribution, with two clear, small maxima, and 
weight also at g ~ 0. 

In the sample of figure (||) instead the e = 1 mean field 
model shows two broad maxima. Again the evolution is 
chaotic, down to the e = 3D model where there are 
four sharp maxima. We want to notice that, if any, the 
structure of the individual Pj{q) becomes more complex 
when approaching the pure 31? model: structures with 
many sharp peaks appear on the larger lattices more in 
the 3-D model than in the mean field one. One reason 
for that is probably connected to the fact that finite size 
effects are larger in the mean field model. Recently a de- 
tailed discussion of models that show a "modified droplet 
model" structure has been reported in [^6). This is very 
interesting since provides a rigorous construction of mod- 
els with a complex structure of states, but looks different 
from what we find here. Pj{q) like the one we observe 
at e = in figure (||) would not appear in a "modified 
droplet model" structure: our results hint that we are 
dealing with a complete RSB like structure of states. 

In figure we show E{q^) as a function of e for L = 8 
and L = 12. The two lines are the best fit of the data to 
a third order polynomial. The most part of this (smooth) 
dependence over e can be related to the renormalization 
of T. The L = 12 data show an even smaller dependence 
over e than the L = 8 data. We have also analyzed the 
fourth and the sixth momentum of the overlap q, that 
have the same smooth behavior. 

We acknowledge valuable conversations with Giorgio 
Parisi. 
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FIG. 1. P^q) as a function of e (going from to 1) and q 3 ^^^^ p.^^^ ^ ^^^^^^ realization 

(from —1 to +1). L = 8, T = 0.7, disorder average taken over „ , , 1 1 i- 1 

^ ' ' ' o 01 trie quenclied disorder. 
150 realizations. 
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